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Abstract. The photochemistry of 32 neutral and 21 ion species in Triton’s atmosphere is 
considered. Parent species N 2 , CH 4 , and CO (with a mixing ratio^of 3 x 1 0 4 in our basic 
model) sublime from the ice with rates of 40, 208, and 0.3 g/cm 2 /b.y., respectively. 
Chemistry below 50 km is driven mostly by photolysis of methane by the solar and 
interstellar medium Lyman-alpha photons, producing hydrocarbons C 2 H 4 , C,H 6 , and C 2 H 2 
which form haze particles with precipitation rates of 135, 28, and 1.3 g/cur/b.y., respec- 
tively. Some processes are discussed which increase the production of HCN (by an order 
of magnitude to a value of 29 g/cnr/b.y.) and involve indirect photolysis of N 2 by 
neutrals. Reanalysis of the measured methane profiles gives an eddy diffusion coefficient 
K = 4XK) 5 cnr/s above the tropopause and a more accurate methane number density near 
the surface, (3.1 ± 0.8)* 10" cm. Chemistry above 200 km is driven by the solar EUV 
radiation (X < 1000 A) and by precipitation of magnetospheric electrons with a total 
energy input of 10 s W (based on thermal balance calculations). The most abundant 
photochemical species are N, H 2 , H, O, and C. They escape with the total rates of 
7.7* 10 24 s \ 4.5* 10 25 s’, 2.4* 10 25 s ’, 4.4* 10 22 s' 1 , and 1.1 *10 24 s' 1 , respectively. Atomic 
species are transported to a region of 50-200 km and drive the chemistry there. Iono- 
spheric chemistry explains the formation of an E region at 150-240 km with HCO as a 
major ion, and of an F region above 240 km with a peak at 320 km and C as a major 
ion. The ionosphere above 500 km consists of almost equal densities of C and N + ions. 
The model profiles agree with the measured atomic nitrogen and electron density profiles. 
A number of other models with varying rate coefficients of some reactions, differing 
properties of the haze particles (chemically passive or active), etc., were developed. 

These models show that there are four basic unknown values which have strong impacts 
on the composition and structure of the atmosphere and ionosphere. These values and 
their plausible ranges are the CO mixing ratio f co = 10 4 -10 3 , the magnetospheric electron 
energy input (1 ± 0.5)* 10 8 W, the rate coefficient of charge-exchange reaction N 2 + + C k 
= lO'MO 10 cm 3 /s, and the ion escape velocity v, 53 150 cm/s. 


1. Introduction 

There are few basic experimental facts about the atmo- 
sphere of Triton. Ground-based spectroscopy shows the 
presence of N 2 , CH 4 , CO, and C0 2 in the surface ice with 
relative fractions 1:5* 10' 4 : 10" 3 : 10 3 [Cruikshank et al , 1984, 
1993; Cruikshank and Apt , 1984], and these species, except 
C0 2 , should exist in the atmosphere. Voyager radio occulta- 
tions [Tyler et al , 1989; Yelle et al. , 1991] measured the 
atmospheric pressure p = 13.5±1 dyn/cm 2 and the mean 
temperature of 40 ± 5 K in the lowest 50 km. The iono- 
sphere with a peak electron density of (3.5 ± 1)* 10 4 cm 3 at 
340 km was observed as well. Voyager Ultraviolet Spectrom- 
eter (UVS) solar occultation spectroscopy [ Broadfoot et al., y 
1989; Krasnopolsky et al, 1993] measured the N 2 density 
profile between 450 and 750 km, which can be described by 
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the N 2 density at 575 km, [N 2 ] = (4 ± 0.4) x 10 8 cm 3 , and the 
isothermal temperature T - 102 ± 3 K. The distribution of 
atomic nitrogen observed by the same instrument corresponds 
to diffusive equilibrium above 300 km with [N] = (1 ± 
0.25)* 10 8 cm 3 at 400 km and T = 100 ± 7 K, and [N] = (5 
± 2.5)* 10 8 cm 3 at 200 km. The total escape rate of atomic 
nitrogen is equal to (1 i 0.3)* 10 25 s 1 [Krasnopolsky et al, 
1993]. Methane and haze profiles have been deduced from 
the Voyager UVS solar occultations [Herbert and Sandel, 
1991]. Both nitrogen and methane are nearly saturated at the 
surface temperature of 38 ± 3 K measured by the Voyager 
1R spectrometer [Conrath et al, 1989] and of 38 ± 1 K 
determined from the profile of the 2.148 jum band of solid 
nitrogen [Tryka et al, 1994]. Haze and clouds seen in 
images of Triton have been studied by Pollack et al [1990] 
and Rages and Pollack [1992]. Study of both UVS and 
television observations of haze yields various properties of 
the haze as functions of height up to 30 km [Krasnopolsky 
et al, 1992; Krasnopolsky , 1993b], of which the specific 
surface is the most important for our problem. Dark plumes 
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(chimneys), which are vertical up to 8 km and which become 
horizontal at this height, indicate the position of the tropo- 
pause [Yelle et al., 1991]. High-energy electrons found in 
Neptune’s magnetosphere in the vicinity of Triton’s orbit 
[Krimigis et al, 1989] are thought to be a source of ioniza- 
tion and heating of Triton's atmosphere. 

These basic facts may be used for more detailed analyses 
of the atmosphere. In spite of some important differences 
between thermal balance calculations by Stevens etal. [1992] 
and Krasnopolsky et al. [1993], both papers showed that the 
mean magnetospheric electron energy input P = 10 8 W was 
required to reproduce the measured densities and tempera- 
tures. Among the differences was an upper limit on the CO 
mole fraction (mixing ratio) f^ < 2x 10 4 given by Stevens et 
al [1992] based on cooling by the CO rotational lines, while 
Krasnopolsky et al. [1993] showed that a CO fraction up to 
10' 2 is compatible with the UVS observations. Later Strobel 
and Summers [1995] came to the conclusion that a model 
with P = 1.4* 10 8 W and 7^ = 2.5 * 1CT 4 provided the best fit 
to the measurements and claimed^ < 3*10' 4 as an im- 
proved upper limit. These revisions were given without 
explanations, and the new limit was too close to their 
preferred value. Besides, they argued that Table 2 in Krasno- 
polsky et al [1993] suggested^ < 10" 3 . However, that table 
was given for fixed values of the magnetospheric electron 
energy input and the heating efficiency. A larger variety of 
models was shown in Table 4 of the same paper, and we do 
not see reasons to reconsider their results. Despite all of the 
differences among the models of Stevens et al 1 [1992], 
Krasnopolsky etal. [1993], and Strobel and Summers [1995], 
the observational constraints made their profiles of tempera- 
ture and density rather close to each other, and the back- 
ground N 2 atmosphere may now be considered as known 
with good accuracy. For example, typical differences 
between the profiles of Krasnopolsky et al. [1993] and of 
Strobel and Summers [1995] are 1-2 K in temperature and 
15-20% in density. 

Few photochemistry papers have appeared since the 
Voyager encounter. Strobel et al. [1990a] considered meth- 
ane chemistry below 100 km, which results in formation of 
the higher hydrocarbons precipitating to the surface, and of 
molecular and atomic hydrogen, which escape. They suggest- 
ed that condensation of the higher hydrocarbons forms the 
haze observed by Voyager. (. Krasnopolsky et al. [1992] use 
this idea as an initial assumption for modeling the haze and 
its UV and visible extinctions and reflectivities.) The first 
ionospheric models were calculated by Majeed et al. [1990] 
and Yung and Lyons [1990], who showed that there is a 
strong coupling between the neutral (lower and upper) 
atmosphere and the ionosphere. Therefore later models by 
Summers and Strobel [1991] and Lyons et al. [1992] treated 
both neutral and ionospheric composition. Lyons et aL [1992] 
pointed out the importance of neutral and ionized atomic 
carbon on Triton. The low ionization potential of carbon 
makes carbon ions stable in the upper atmosphere, which 
mostly consists of molecular and atomic nitrogen and 
hydrogen. Solar radiation is the only source of ionization in 
the model of Lyons et al. [1992], and a very long lifetime of 
C + permits a fit to the ionospheric observations. A charge- 
exchange reaction between N 2 + and C was a source of C + , 
and the required rate coefficient was very large, 10 9 cm 3 /s. 

The discovery of CO in the surface ice [Cruikshank et al. , 
1993] stimulated further studies of Triton’s photochemistry. 


Three preprints appeared in 1993 on this subject: J. R. Lyons 
et al. (Photochemistry of the atmosphere and ionosphere of 
triton, submitted to Icarus , 1993; hereafter referred to as 
LY A93), V . A. Krasnopolsky and D. P. Cruikshank (unpub- 
lished data, 1993), and Strobel and Summers [1995] (hereaf- 
ter referred to as SS95. The presence of CO makes the 
photochemistry more complex and interesting due to the 
appearance of oxygen species. Furthermore, carbon monox- 
ide, like methane, is a source of carbon atoms, and it makes 
it possible to model the ionosphere profile with a more 
realistic rate coefficient of the N 2 + + C reaction. The condi- 
tions of saturation of CO on Triton correspond to the CO 
mixing ratio of 7%. Then the application of Raoult’s law to 
the CO mixing ratio in the ice, la 3 , yields its mixing ratio of 
0.07x 10 3 = 7* 10' 5 in the atmosphere. However, the atmo- 
sphere-ice system may not be in local equilibrium, as is the 
case for methane, for which the atmospheric abundance 
exceeds that predicted by Raoult’s law by 3 orders of 
magnitude. That is why measurements of the CO abundance 
in the atmosphere are badly needed. 

LYA93 will be a basic reference in our study here. They 
assumed the CO mixing ratio of 10 r 4 in the atmosphere and 
succeeded in reproducing the measured peak electron density 
by the solar model (without precipitation of magnetospheric 
electrons) and with the N 2 + + C reaction rate coefficient (k^ 
in our notation) of 10 10 cm 3 /s. SS95 is a review chapter 
which does not discuss many important details of their model 
( e * profiles of only five species from their latest model 
are shown, and other profiles given in the paper represent the 
corrected version of Summers and Strobel [1991]). Their 
model is based on the magnetospheric electron energy input 
of 1 .4* 10 8 W, the CO mixing ratio f co = 2.5 * 1CT 4 , and k% 4 = 
10 11 cm 3 /s. We will compare our models with SS95 where 
possible. 

As with LYA93 and SS95, the objective of this paper is 
a self-consistent model of the atmosphere of Triton including 
its ionosphere. However, our models differ in some important 
details from those papers, and this difference is what has 
stimulated this work. 

2. General Ideas of Photochemistry of the 
Atmosphere and Ionosphere of Triton 

General ideas of the photochemistry of the atmosphere and 
ionosphere of Triton were discussed in previous papers, and 
we descibe them very briefly to give more attention to 
specific problems. Solar EUV radiation shortward of 1000 A 
and precipitating magnetospheric electrons ionize and 
dissociate N 2 and produce N, N 2 \ and N* with an N//N* 
production ratio of 5 (this ratio is almost the same for EUV 
photon and low-energy electron impacts). This occurs mostly 
at 200-500 km. Solar^ radiation in the range 1000-1400 A 
(Lyman-alpha 1216-A photons constitute 80% of this 
radiation) and the Lyman-alpha interstellar background with 
an intensity of 340 Rayleigh (SS95) dissociate methane in 
the lowest 40 km to produce CH 3 , CH, H, and H2. Some part 
of the methane is rebuilt from CH 3 and H, but further 
destruction of CH 4 in a reaction with CH results in a total 
yield of the methane photolysis y « 1 (0.6 in Strobel et al. 
[1990a] and LI in our models). Reactions of CH 3 and CH 
form more complex hydrocarbons; C 2 H 4 , C 2 H 6 , and C 2 H 2 are 
the most abundant. The saturation densities of these hydro- 
carbons are of the order of a few molecules per cubic 
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centimeter at 40 K, and therefore these species condense and 
form a haze, which was observed up to 30 km by the 
Voyager 2 TV and UVS instruments. The ratio C:H = 1 :4 in 
CH 4 and * 1:2 in the photochemical ly formed hydrocarbons, 
and therefore 2y hydrogen atoms are released per photon due 
to the hydrocarbon chemistry. These hydrogen atoms 
constitute an upward hydrogen flow, which escapes and 
which is equal to the 1000-1400 A photon flux multiplied by 
2y, i.e. 5x10 s cm'V 1 and 1 . 1 5 x 1 0 26 s' 1 from the whole 
satellite. The H 2 :H ratio is determined in this flow by 
chemical processes in the atmosphere and ionosphere, 
however, the total hydrogen flow is rather stable and should 
not depend on a chosen chemical model. This flow is 
diffusion limited, so the (2^ + H) mixing ratio should be 
rather constant with height in the atmosphere and equal to / 
» <&H/b * 2* 10" 4 . Here is the flow, H * 15 km is the scale 
height, b = Dn = LSSxlO 17 ^ 82 cm'V 1 [ Chamberlain and 
Hunten, 1987], and D is the diffusion coefficient of H 2 in N 2 . 
The loss of methane by photolysis is compensated by its 
transport by eddy and molecular diffusion, and the measured 
methane profiles [Herbert and Sandel, 1991] are used to 
determine the eddy diffusion coefficient. 

The neutral atmosphere at ionospheric heights, i.e., above 
200 km, consists of molecular and atomic nitrogen and 
hydrogen and CO, which is rather inert below 200 km. The 
primary N 2 + and N + ions react with H 2 to form N 2 H + and 
NH* and recombine. Some part of the primary ions exchange 
their charge with CO and form CO\ Recombination of CO + 
is the main source of C and O atoms in the atmosphere. 
Small fractions of the N, C, and O atoms escape, and 
downward flows of these atoms interact with each other and 
with other components of the atmosphere below 150 km. The 
primary result of this interaction is recombination to N 2 and 
CO, though some HCN, other nitriles, H 2 CO, and C0 2 form 
and condense. 

The lifetime of N 2 + is rather short due charge-exchange 
and recombination processes, and its density should be far 
smaller than the measured peak electron density (3.5 ± 
1)* 10 4 cm' 3 . Radiative recombination of N + is much 
slower than dissociative recombination of N 2 + and other 
molecular ions. However, FT can be lost in the charge- 
exchange reactions with H 2 and CO. Therefore the first 
ionospheric models [Majeed et ai , 1990; Yung and Lyons , 
1990] explored an idea that the ion column production rate 
could exceed the H 2 production in the lower atmosphere, and 
the excess IT ions should recombine radiatively, forming a 
dense ionosphere. Later, Summers and Strobel [1991] 
reduced the requirement to the ion production rate by 
enhanced formation of H + in charge-exchange of N 2 + with H. 
W vanished in reactions with hydrocarbon and nitriles. Later 
measurements gave only a rather low upper limit to a rate 
coefficient of this reaction, and H + densities are strongly 
reduced in the case of CO (see below). In the presence of 
CO and hence C, the charge-exchange reaction between N 2 + 
and C forms C + ions. These ions have a very long chemical 
lifetime because they can either recombine radiatively or 
react with HCN, CH 4 , and the other hydrocarbons, which are 
very scarce at the ionospheric heights. 

These are the basic ideas of Triton’s photochemistry. The 
experimental N 2 densities and thermal balance calculations 
provide the background N 2 atmosphere. The measured CH 4 
densities are used as the boundary condition for this parent 
species and allow us to determine the eddy diffusion coeffi- 


cient. The measured electron density profile and the atomic 
nitrogen profile are used to check a validity of a model. 

3. Reactions and Their Rate Coefficients 

One of the problems in modeling Triton’s photochemistry 
is the lack of measurements of reaction rate coefficients at 
temperatures typical of Triton's atmosphere. The situation is 
fairly good for ionospheric reactions which typically exhibit 
small temperature dependences. The temperature behaviors 
have been measured for some basic ion reactions (see 
references in work by Anicich [1993]). Our knowledge is less 
satisfactoiy for neutral reactions in the upper atmosphere 
where T = 100 K, and is very poor for reactions in the lower 
atmosphere with T * 40-50 K. 

Some reactions proceed via an energy barrier E (in 
Kelvins), and their rate coefficients are given in the Arrhe- 
nius form k = A e' m . Temperature behavior of a reaction 
without the energy barrier is determined mostly by a poten- 
tial energy function of interacting molecules and may be 
given as k = AfiOQITJF. Both types may combine to form k 
= A(3QQIT) a & E/r i which is used, for example, in the RRKM 
approximation. Our reference guide for neutral reactions is 
Chemical Kinetics Database [Mallard et ai, 1994]. The 
detailed physics is unknown in many experimental studies of 
reaction rate coefficients, and one of the above functions is 
sometimes arbitrarily chosen to represent the experimental 
results. Typically this choice is not important in a measured 
temperature range, but it is crucial for extrapolation to 
Triton’s temperatures. Therefore we transform rate coeffi- 
cients of all three-body association reactions, given in the 
Arrhenius form, to the power form, because these reactions 
proceed typically without an energy barrier. For example, a 
rate coefficient of 

O + O + N 2 — * 0 2 + N 2 

was measured by Campbell and Gray [1973] and given as k 
= 9.5* lO-V 85 * cm 6 /s. This corresponds to k = 5*10“* 
(300/7) 2 cm 6 /s in the temperature range of their measure- 
ments, 196-300 K, and this form is widely used in photo- 
chemistry. Furthermore, we are cautious in extrapolations of 
rate coefficients given with a negative activation energy E , 
and either the power form or constant values will be used. 

The number of compounds consisting of N, H, C, and O 
atoms is so huge that we try to reduce both the number of 
compounds and their reactions in our model to a reasonable 
minimum. Shown in Table 1 are 75 reactions of 31 neutral 
species, and in Table 2, 63 reactions of 20 ions. We show in 
these tables the main reactions which determine balances of 
the parent species (N 2 , CH 4 , CO), the principal photochemi- 
cal products (N, H, H 2 , C, O, C 2 H 4 , HCN, C 2 H^, and 

the principal ions (C + , N\ N 2 \ N 2 H + , CO + , and HCO + ). 
Though hydrocarbons and HCN are very scarce at ionospher- 
ic heights, we take special care to describe in detail their 
reactions with ions because their densities, even of the order 
of 1000 cm' 3 , may be important for C\ To minimize the 
number of species and their reactions, we do not use reac- 
tions which do not change production and loss of the basic 
species. For example, the net effect of 

HCCT + HCN - HCNH + + CO 
HCN1T + e -» HCN + H 

net HCO + + e - CO + H 

adds almost nothing to dissociative recombination of HCO\ 
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Table 1 . Neutral Reactions in Triton’s Atmosphere, Their Rate Coefficients and Column Rates 



Reaction 

Rate Coefficient 

Reference 

Column Rate 

la 

lb 

CH 4 + hv -* CH, + H 
-*■ CH + H + Ho 

y^O.49 7.9*10" 
0.506 

Mordaunt et al. 1 1993] 

1.10+8 T 

1.14+8 

lc 

C + H. + H, 

0.004 

Braun et al. [ 1 969 1 

9.0045 

2 

CTE + hv - C.H + H 

9*10'' 

Okabe |1983| 

2.2444 

3 

C,H 4 + hv - CoH, + H + 0.5H, 

3* Iff® 

Yung et al. (1984] 

5.6045 

4a 

C 7 H 6 + hv - C.H 4 + H + 0.5 Ho 

4* 10' 9 

Yung et al. 1 1984] 

1.3544 

4b 

- CoHo + 2H. 

2* 10' 9 

Yung et al. [1984] 

6.7743 

5 

HCN + hv - H + CN 

2 x 1 0 8 

Yung et al. (1984] 

2.8045 

6 

HCNo + hv CH + No 

3* Iff 7 

see text 

1.2646 

7 

N, + hv N + N 

- 

Zipf and McLaughlin ( 1 978] 

1.77+7 

8 

CO + hv — C + O 

- 

see text 

3.80+5 

9 

N 2 + £-»N + N + tf 

- 

Fox and Victor [1988] 

5.04+7 

10 

CO + e->C + Q + e 

- 

similar to No 

1.64+3 

11 

N + N+ M-*No + M 

4* 10'”(300/7) 2 

Campbell and Thrush [1967] 

1.67+7 

12 

H + H + M— >Ho + M 

9* 10' "(300/7)' 3 
5.5* 10' 31 (300/7)' 6 

Tsang and Hampson [ 1986] 

6.85+7 

13 

C + C + M-*-Co + M 

Slack 1 19761 

6.23+4 

14 

O + O + M— ►Oo + M 

5 * 10'”(300/7) 2 

Campbell and Gray [1973] 

1.54+5 

15 

c + o + m->c6 + m 

5*10"(300/7) 2 

assumed 

7.56+3 

16 

N + H + M-NH + M 

5*10 ,2 (300/7) 

Brown [1973J 

3.60+7 

17 

N + C + M CN + M 

10 ,2 (300/7) 

Kley et <?/.[! 974] 

3.42+4 

18 

N + O + M NO + M 

lff ,2 (300/7) ofis 

Campbell and Gray [1973] 

2.87+5 

19 

H + O + M — ► OH + M 

4.4* 10 ,2 (300/7) 

Tsang and Hampson [1986] 

2.04+5 

20 

H + CH, + M - CH 4 + M 

6*10' 29 (300/7)‘ 8 

Baulch et al. [1992] 

7.18+7 

21 

Ho + C + M - CHo + M 

1.4*10 31 (300/7) 

Husain and Young [1975] 

2.91+5 

22 

Ho + CH + M - CH, + M 

6* lO' ,0 (3OO/7) 2 3 

Berman and Lin [ 1984] 

4.20+5 

23 

CH, + CH, + M — CoH 6 + M 

10 23 

La ufer et al. [1983] 

1.78+7 

24 

CH + No + M -> HCNo + M 

2.6*10’ 31 (300/7) 

Berman and Un [1983a] 

2.14+7 

25 

C + No + M - CNN + M 

3 *10 ”(300/7) 

Husain and Kirsch [1971] 

5.05+7 

26 

CH 4 + CH -+ CH 4 + H 

5* 10' n e 200/r ; 10‘ 10 * 

Berman and Un [1983b] 

9.46+7 

27a 

CH, + N -*■ HoCN + H 

4.3*10‘ l V 42a " 

Mars ton et al. [1989] 

3.15+5 

27b 

- HCN + Ho 

4«lV n e 420 ' T 

Mars ton et al. [ 1989] 

2.94+4 

28 

CH, + O HoCO + H 

1.4*10 10 

Baulch et al. [ 1 992] 

3.19+6 

29 

CHo + H - CH + Ho 

10‘ 10 * 

Baulch et al. [1992] 

3.33+6 

30 

CHo + N -> HCN 4- H 

5*10 n e 250fT 

Yung et al. [1984] 

4.34+4 

31 

CHo + O - CO + Ho 

2* 10 10 

Baulch et al. [ 1 992] 

1.43+6 

32 

CH + N-> CN + H 

2* 10 11 

Messing et al. [1981] 

1.91+6 

33 

CH 4 H - C 4 Ho 

5 *10 11 

Harding et al. [1993] 

2.82+6 

34 

CH 4 C - Co 4 H 

5* Iff" 

assumed 

2.05+4 

35 

CH + O -* CO + H 

5*10“ 

Baulch et al [1992] 

3.27+5 

36 

CH 4 HoCO - CHo 4 HCO 

10 10 * 

Baulch et al. [1992] 

2.99+6 

37 

NH 4 H -* N 4 Ho 

5*10“ 

Morle y [1981] 

4.58+7 

38 

NH + N -► No + H 

5*10" 

Miller et al [1985] 

1.03+7 

39 

NH 4 C - CN 4 H 

5* 10 11 

assumed 

2.37+3 

40 

NH 4 0*-+ NO 4 H 

1.2* 10 10 

Cohen and Westberg [1991] 

4.29+5 

41 

CNN 4 N - CN 4 N 2 

2*10" 

assumed 

4.13+7 

42 

CNN 4 H ^ CH 4 No 

2*10" 

assumed 

5.64+6 

43 

CNN 4 C * Co + No ~ 

2*10" 

assumed 

2.69+5 

44 

CNN 4 O - CO + No 

2*10" 

assumed 

3.35+6 

45 

CN 4 N -> C 4 No 

10 10 

Whyte and Phillips [1983] 

4.53+7 

46 

CN + O -> CO + N 

1.7*10" 

Baulch et al [1992] 

6.16+5 

47 

CN 4 CoHo -> HCN 4 CoH 

2.7*10 lo (300/7)°V 19 ' r 

Sims et al [ 1993] 

6.22+4 

48 

CN 4 CoH 4 - HCN 4 CH, 

2.5* 10 lo (77300)°V 28,r 

Sims et al [1993] 

2.94+4 

49 

CN 4 CH 6 - HCN 4 CH, 

2.1*10"(r/300) 022 / 8 " 

Sims et al. [ 1993] 

3.39+3 

50 

CN 4 HoCO - HCN 4 HCO 

5*lff l4 7* 

Yu et al [1993] 

2.07+5 

51 

CoH 2 + C — * C, + Ho 

10"* 

Haider and Husain [1993] 

1.02+6 

52 

CoH 4 4 c - CoHo + CHo 

10"* 

Haider and Husain [1993] 

1.47+6 

53 

C4I 6 4 C * CH“ 4 CHo 

Iff" 

assumed 

5.80+4 

54 

C 2 H 4 4 O — * CH, 4 HCO 

1.5*10' l2 (77300) 1 V 2l5/r 

Tsang and Hampson [1986] 

1.70+5 

55 

N 4 o - NO 4 //v 

1.8*10' l7 (300/7)° 5 

Dalgartio et al [1992] 

4.96+7 

56 

N 4 oh -4 NO 4 H 

5*10"(300/7) 025 

Howard and Smith [1981] 

1.99+5 

57 

N 4 NO - No 4 O 

3*10" 

Baulch et al [1992] 

5.06+7 

58 

CO + OH -► COo 4 H 

2.6*10 ‘V wr 

Frost et al [1993] 

5.30+3 

59 

HoCN 4 N HCN 4 NH 

iff'V 200 ' 1 ' 

Marston et al [1989] 

3.15+5 

60 

HCNo 4 H - HCN 4 NH 

10' 4 

assumed 

2.02+7 

61 

Oo + C — ► CO + O 

10' 12 * 

Dor the et al [1991] 

1.54+5 

62 

C 2 4 o - CO 4 c 

10 12 * 

Fairbairn [1969] 

9.52+4 

63 

Co 4 N -4 CN 4 C 

10 12 

assumed 

1.27+6 

64 

C, 4 O -4 CO 4 Co 

10 12 

assumed 

7.09+4 

65 

C, 4 N -> CN 4 c; 

10 12 

assumed 

9.46+5 

66 

C 2 H 5 + H -* CH, + CH, 

6*10" 

Baulch et al (1992] 

3.39+3 
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Table !• (continued) 



Reaction 

Rate Coefficient 

Reference 

Column Rate 

67 

CoH 3 + H - C.H, + Ho 

2* 10“ 

Baulch et al. 1 1992] 

2.94+4 

68 

C.H + O - CO + CH 

1.7*10-" 

Baulch et al. [ 1 992] 

3.97+3 

69 

CJf + N — CN + CH 

10" 

assumed 

8.05+4 

70 

C^H + CH 4 -» C,H, + CH, 

3*10 l V 25<vr 

Tsang and Hampson (19861 

134 

71 

HCO + H — CO + H, 

1.5* 10 10 

Baulch et al. ( 1 992] 

3.19+6 


Photolysis rate coefficients (in s 1 ) are halves of those on Neptune’s orbit for X > 80 nm. Rate coefficients are in cmVs 
and cm 6 /s for second- and third-order reactions, respectively, and column rates (in cm 2 s l ) are normalized to the surface. 
* Values adopted for Triton's conditions. 
f 1.10+8 = l.lOxlO 8 . 


The same holds for some branches of dissociative recombina- 
tions and ion-molecule reactions, e.g., 

n 2 + + o - cr + n 2 

Q + + N, — NO* + N 
net N 2 + + O - NO* + N 

just repeats the main branch of the N 2 + + O reaction. Before 
we neglect any species or reaction, we check qualitively for 
possible effects of this neglection. 

Now we begin to discuss some reactions and their rate 
coefficients which are new or differ substantially from the 
previous models. 

Photolysis of methane was measured recently by Morda- 
unt et al [1993] to proceed via 

(Rl) CH 4 + hv - CH 3 + H y = 0.49 

CH + H + H 2 y = 0.5 1 

Previously it was thought that CH 2 is the main photolysis 
product (y = 0.92), the rest being CH [e.g., Atreya 1986], 
Then reactions 

CH 2 + H -+ CH + H 2 
CH + CH 4 - C 2 H 4 + H 

increased the effective photolysis yield to 1.75, and forma- 
tion of ethane and acethylene was extremely low [Krasnopol- 
sky etaL, 1993]. The high yield of CH 3 increases three-body 
formation of CH 4 (reducing the effective photolysis yield to 
1. 1) and ethane QH^ 

It was argued by Krasnopolsky et. al [1993] that a rate 
coefficient of 

(R29) CH 2 + H -* CH + H 2 

given as £29 = 4.7* 10* 10 exp(-370/7) cm 3 /s (e.g., LYA93) is 
not applicable to the low Triton temperatures (it is equal to 
5* 10 14 cm 3 /s at 40 K). Baulch et al.' [1992] recommend k 29 
= 10 10 exp(300/7), which is too big at low temperatures, and 
a value of l O' 10 cm 3 /s will be used. This reaction was of great 
importance for the old scheme of the methane photolysis 
which produced much CH 2 . Its impact is much smaller now. 

Current photochemical models presume a formation of 
HCN on Triton by (R30) N+CH 2 and (R27) N+CH 3 and 
ignore the fact that CH radicals react quickly with N 2 , the 
reaction having been observed many times in laboratory (19 
references to measurements of this reaction are given by 
Mallard et al. [1994]). It is thought that this reaction is 
responsible for dissociating N 2 molecules in flames and 
combustion. However, the reaction to form HCN, 

CH + N 2 — HCN + N 

is slightly endothermic and cannot proceed even at room 
temperature. The most detailed analysis of this process was 
made by Berman and Lin [1983a], who studied it over a 


large range of temperatures and pressures. The very complex 
behavior of the system suggests that the initial process at 
room and low temperatures is 

(R24) CH + N 2 + M - HCN 2 + M k 24 = 2.8* 10' 31 cm 6 /s 

This reaction may be followed by a reaction of HCN 2 , with 
CH forming two HCN molecules with a net energy release 
of 9.44 eV. N and H atoms are much more abundant on 
Triton than CH, and these are possible candidates to react 
with HCN 2 . The rate coefficient k 24 is rather high; therefore 
(R24) is spin-allowed and HCN 2 has spin 1/2. Then the 
reaction with N( 4 S) to form HCN+N 2 is spin forbidden, and 
(R59) H+HCN 2 is the most probable reaction to proceed. 
Products other than HCN+NH may form, but we ignore other 
branches and adopt k m = 10' 14 cm 3 /s for this branch. Berman 
and Lin mentioned that the HC-N 2 bond energy is 2 eV and 
HCN 2 dissociates by visible light with a cross section close 
to that of diazomethane. The latter would have a dissociation 
rate coefficient of 3*10" 7 s' 1 on Triton (cross section from 
Okabe [1978]). Reactions (R24) and (R60) are the main 
processes of HCN formation in our model. Berman and Lin 
[1983a] predict an increase of L, 4 by more than an order of 
magnitude from 300 K to 100 K, though they do not 
preclude the possibility of a small activation energy. The 
value chosen here reflects these considerations. 

The formation of NH in (R60) is mostly followed by 
(R37) NH + H. Therefore (R24, 60, and 37) may be consid- 
ered as a chemical cycle: 

CH + N 2 + M -*• HCN 2 + M 

HCN 2 + H - HCN + NH 

NH + H -> N + K, 

net CH + N 2 + 2H - HCN + N + H 2 

which forms HCN, accelerates the association of H atoms to 
form H 2 , and provides a rare case of indirect photolysis of N 2 
by Lyman-alpha radiation (which produces CH from CH 4 ). 
These species are among the most abundant photochemical 
products. To check the effect of this cycle, we will consider 
a model without this cycle and compare it with a similar 
basic model. 

Transition probabilities of two components of the N(*D 
~* 4 S) 5200 A doublet are 

N( 2 £>j^ - N( 4 S) + hv A m = 1.6*ia 5 s' 

N( 2 Dj „) - N( 4 S) + hv A m = 0.7 x ia s s' 1 

[e.g., Allen, 1973], Then a weighted-mean transition proba- 
bility of N^D) is equal to 1.06*1 0' 5 s' 1 using statistical 
weights g = 2J + 1. We discuss this to prevent a further 
spread of an erroneous value 2.3* 10" 5 = 1.6* 10 5 + 0.7* 10" 5 
in photochemical papers (e.g., Stevens et al. [1992], LYA93, 
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Table 2. Reactions in Triton’s Ionosphere, Their Rate 
Coefficients, and Column Rates Normalized to the Surface 


Reaction Rate Coefficient, Column Rate, 

cm 3 /s cm 2 s 1 


72a 

N, + - N 2 * + 2e 

- 

5.64+T 

72b 

- N* + N + 2e 

- 

1.24+7 

73a 

CO + e -* CO* + 2e 

- 

2295 

73b 

— C* + O + 2e 

- 

503 

74 

N + e — N* + le 

- 

6.84+5 

75 

C + e — C + le 

- 

1.76+4 

76 

O + e -* 0 + + 2e 

- 

2.35+4 

77a 

No + /?v — No* + e 

- 

2.92+7 

77b 

- N* + N + e 

- 

4.17+6 

78a 

CO + hv -*■ CO* + e 

- 

739 

78b 

-*■ C* + O + e 

- 

116 

79 

N + hv - N* + <? 

- 

3.88+5 

80 

C + hv - C* + e 

- 

6.91+4 

81 

O + hv —*■ 0 + + e 

- 

1.30+4 

82 

N 2 * + Ho - N,H* + H 

1.3 xlO' 9 * 

2.25+7 

83 

n 2 * + co -*• co + + n. 

7.3*10-“ 

4.78+5 

84 

N 2 * + C - C* + No 

4*10 11 * 

1.64+6 

85 

N 2 * + O - NO* + N 

1.4* 10 10 

9.33+6 

86 

N 2 * + N - N* + N 2 

10" 

1.75+7 

87 

n 2 * +ch 4 - ch 3 * + h + n 2 

1.2 *10‘ 9 * 

1.67+4 

88 

N 2 + + HCN - HCN* + N 2 

3.9* 10 10 

586 

89 

N* + H 2 - NH* + H 

3*10 10 * 

1.94+7 

90 

N* + CO -*• CO* + N 

9*10 10 * 

1.48+7 

91 

N* + CH 4 - CH 3 * + NH 

10 9 * 

9.34+3 

92 

N* + HCN — HCN* + N 

3.7* 10- 9 

4.45+3 

93 

CO* + H 2 -+ HCO* + H 

1.4* 10‘ 9 

4.27+6 

94 

CO* + H -+ H* + CO 

7.5*10 ! ° 

3.86+5 

95 

CO* + C - C* + CO 

4*10 " * 

1.90+5 

96 

CO* + O - O* + CO 

1.4*10 10 

1.21+6 

97 

CO* + ch 4 - ch 4 * + CO 

1.34* 10' 9 

4.91+3 

98 

CO* + HCN -+ HCN* + CO 

3.4 *10' 9 

1.36+3 

99 

H* + O - O* + H 

8.5*10"* 

3.92+5 

100 

H* + CH 4 - CH 3 * + H 2 

4. 15 *10 9 

8.01+3 

101 

H* + HCN - HCN* + H 

1.1*10^ 

2.47+3 

102 

O* + N 2 - NO* + N 

2.8 *10 12 * 

1.53+6 

103 

O* + H 2 - OH* + H 

1.6* 10 9 

9.75+4 

104 

O* + H - H* + O 

6* 10 10 * 

7.71+3 

105a 

C* + CH 4 - cyv + H 

9.4*10' 10 

2.72+5 

105b 

- C.H.* + H, 

3.6* 10 10 

1.03+5 

106 

C* + HCN -+ C, N* + H 

3*10' 9 

1.76+5 

107 

NH* + No - N 2 H* + N 

6.5*10 10 

1.94+7 

108 

N 2 H* + CO - HCO* + No 

8.8 *10'° 

5.61+6 

109 

N 2 H* + CH 4 - CH 5 * + Nz 

8.9 *10 10 

2.62+4 

110a 

C,H 2 * + N - CHCN* + H 

1.5*10 10 

6.23+4 

110b 

-> C,N* + H2 

7.5*10-“ 

3.12+4 

110c 

— CH* + HCN 

2.5*10" 

1.04+4 

111 

Cft* + N — CHCN* + H 2 

2.2*10 10 

2.72+5 

112 

OH* + N 2 - N 2 H* + O 

2.4* 10 10 

9.75+4 

113 

C* + e -+ C + hv 

8*10 12 * 

2.26+5 

114 

N* + e - N + hv 

7* 10 12 * 

4.02+4 

115 

N 2 * + e - N + N 

5.7*10-®* 

3,41+7 

116 

CO* + * - c + o 

1.7* 10‘ 7 * 

9.19+6 

117 

NO* + e -> N + O 

6.5* 10 7 * 

1.09+7 

118 

N 2 H* + * - No + H 

6.3 *10' 7 * 

3.64+7 

119 

HCO* + e — H + CO 

5*10" 7 * 

9.89+6 

120 

HCN* + <> — H + CN 

3.5 * 10" 7 * 

8.86+3 

121a 

CHCN* + e -* CH + CN 

2.6* 10' 7 * 

1.67+5 

121b 

— C + HCN 

2.6* 10 7 * 

1.67+5 

122 

cysr + e -+ CN + C 

5*10 7 * 

2.07+5 

123 

CH* + e -> C + H 

5*10" 7 * 

1.04+4 

124 

CH 3 * + e - CHo + H 

5*10" 7 * 

3.41+4 

125 

CH 4 * + <> - CH 3 + H 

5*10" 7 * 

4.91+3 

126 

CH 5 * + * - CH 3 + Ho 

5*10" 7 * 

2.62+4 


Rate coefficients for ion-molecular reactions are taken from 
Anicich [1993], and their temperature dependences (if measured) 
from references therein. Radiative recombination coefficients are 
from Prasad and Huntress [1980], and dissociative recombination 
coefficients are from Mitchell [1990] with some values assumed. 

* For T = 100 K. 

f 5.64+7 = 5.64* 10 7 . 

* The value assumed for the basic model. 


SS95). Because chemical quenching of N( 2 D) is small or 
nonexistent under the conditions on Triton, we ignore 
chemical effects of N( 2 /)). 

The formation of CNN (R25) and its reaction with N is 
the important chemical cycle of atomic nitrogen loss 
(LYA93, SS95): 

(R25) C + N 2 + M -+ CNN + M 
(R41) CNN + N - CN + N 2 
(R45) CN + N -» C + N, 
net N + N -> N 2 

LYA93 assumed k M and rate coefficients of CNN with H, O, 
and C as 10 ~ 10 cm 3 /s based on comparison with (R45), while 
SS95 adopted 10' n cm 3 /s. We consider CCO as the best 
analog for CNN due to the similarity of the CO and N 2 
molecules. The C-CO bond energy is 2.2 eV [ Okabe , 1978], 
and reactions of CCO with O and H are really in the range 
of 10 n -10 10 cm 3 /s [Bauer et al, 1985]. We will use k = 
2* 10 ' 11 cm 3 /s for the reactions of CNN with N, H, C, and O 
and ignore the possible interactions of CNN with hydrocar- 
bons adopted in LYA93. Note that if the C-N 2 bond energy 
exceeds 1.6 eV, then the reaction branch 

CNN + H - CN + NH 

is endothermic. Therefore we ignore this branch (suggested 
by LYA93 and SS95) and assume CH + N 2 as the only prod- 
ucts of the CNN + H reaction. N-N bonds are very strong; 
therefore, breaking of these bonds may be a minor branch in 
the CNN + C reaction. Therefore we suggest 
(R43) CNN + C-C 2 + N 2 

instead of CNN + C — * ► CN + CN suggested by LYA93. 

Reactions (R47-49) of CN with hydrocarbons [Sims et al , 
1993] were observed recently at temperatures down to 25 K. 
The branching ratio of C 2 H 3 CN is equal to 0.2 in the reaction 
CN + C 2 H 4 [Monks et al, 1993], while that of HC 3 N in 
(R47) CN + C 2 H 2 is considered negligible by Sayah et al 
[1988]. HCN is the only nitrile considered in all of the 
experimental studies of CN + C 2 H 6 known to us. Therefore 
we assume HCN is the only nitrile in these reactions and 
neglect the possible production of others. 

Recent studies of reactions (R51) and (R52) of C with 
C 2 H 2 and C 2 Hi by Haider and Husain [1993] are important 
for Triton’s chemistry. The rate coefficients were measured 
at 300 K, and an order-of-magnitude reduction is adopted by 
us to apply them to Triton’s conditions. The products were 
not identified in these reactions, and we suggest the most 
exothermic branches. 

A rate coefficient of the key reaction in the C + ionosphere 
(R84) N 2 + + C - C + + N 2 + 4.32 eV 

is unknown. Based on its similarity to 

N 2 + + H - H + + N 2 + 2.0 eV, 

SS95 suggested k% 4 » 10 11 cm 3 /s. We agree that k% A may be 
as low as 10 ‘ 11 , but do not see good reasons to preclude 
higher values. Consider, e.g., 

(R96) CO + + O -► 0 + + CO + 0.4 eV k^ = 1.4* 10 10 cm 3 /s 
[Fehsenfeld and Ferguson y 1972; Anicich , 1993]. The oxygen 
atom is a much better analog to C in terms of thermal 
velocity and electron state than the hydrogen atom is, and 
CO + is a good analog to N 2 + . The high exothermicity of 
(R84) does not preclude a high rate coefficient, e.g., 

H + + NO — NO + + H + 4.35 eV k= 1.9* 10 9 cm 3 /s 
[Fehsenfeld and Ferguson , 1972; Anicich , 1993], We assume 
i ^4 = 4*10‘ n cm 3 /s in our basic model. The same value is 
adopted to a reaction 
(R95) CO* + C C* + CO 
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Some models with smaller and larger values of k% 4 and k 95 
will be considered as well. 

Prasad and Huntress [1980] recommend values of 
4.4* 10* 12 (300/7)° 6 and 3.8* 10 12 (300/7)°* cmVs for radiative 
recombination coefficients of C + and N + , respectively. We 
adopt these values, which are larger by a factor of 1 .5 at T 
= 100 K than 5* 10' 12 cmVs adopted by LYA93. The smaller 
radiative recombination coefficients facilitate a fit to the 
measured electron densities. 

4. Model 

Background atmosphere and energy sources. We use 
in our model the profiles of temperature and N 2 density from 
Krasnopolsky et ai [1993]. Our model is steady state and 
one-dimensional, and reflects global-mean and orbit-averaged 
conditions at high solar activity typical of the period before 
the Voyager 2 observations (F l01cm « 200). Three sources of 
dissociation and ionization are considered: the solar EUV 
radiation, the interstellar background Ly man-alpha radiation 
with a mean intensity of 340 R (SS95), and precipitation of 
magnetospheric electrons with a total energy input of 10 8 W. 
Solar EUV intensities are taken from Torr and Torr [1985] 
and interpolated to the solar activity index F l 07cm = 200. We 
scale the H Lyman a flux 4. 1 * 10 8 photons cm V 1 from Lean 
[1991] by a factor of 1.28 to find the flux in the 1000-1400 
A interval. The magnetospheric electron ionization profile is 
taken from Strobe l et ai [1990b], is shifted upward by two 
scale heights according to Summers and Strobel [1991], and 
is multiplied by a factor of 0.162, which reflects a ratio of 
the global-mean and orbit-averaged electron flux to that near 
Neptune’s magnetic equator in the chosen model of Krasno- 
polsky et al [1993]. The photolysis rates are halves of those 
calculated for solar zenith angle of 60°. The CO mixing ratio 
is 3 x 10" 4 in the basic model. This value is just in the middle 
of the 10 M0 3 range, which may be the most plausible. 
SS95 prefer 2.5* 10" 4 , which slightly exceeds their previous 
upper limit < 2* 10' 4 [Stevens et ai , 1992] and is very 
close to their current upper limit f co <3* lO^ 4 . Uncertainties 
in the magnetospheric electron energy input and heating 
efficiencies of some processes are rather large, and we do 
not understand why the difference between the preferred CO 
mixing ratio and its upper limit is so small, with both values 
based on the thermal balance calculations. 

Haze model. An important feature of chemical modeling 
of Triton is the condensation of hydrocarbons (except CH 4 ) 
and nitriles on the haze particles. Saturation densities of 
these species are extremely low at T « 40 K (~10 cm 3 for 
C 2 H 4 [Brown and Ziegler ; 1979]), and each collision with a 
particle results in condensation. The condensation rate 
coefficient is therefore 

x c 1 - V^/4 

where v ( is thermal velocity of a molecule, S is the haze 
specific surface, i.e., a sum of surfaces of all particles per 
cubic centimeter. The properties of the haze as a function of 
height were determined by Krasnopolsky et ai [1992] and 
Krasnopolsky [1993b] up to 30 km based on the Voyager 2 
UVS and TV observations [Herbert and Sandel, 1991; 
Pollack et ai, 1990]. From those data, x c _i = 9*10’ 5 exp(- 
h/H) s' 1 . Here H = 14 km is the atmospheric scale height. 
This formula can be used below 30 km, and two scenarios 
are possible above this height. In the first, nucleation and 
spontaneous condensation occur slightly above 30 km, say, 
at 40 km, and no particles exist above this height. The 


second case corresponds to condensation of all locally 
produced material to grow the particles. We consider this 
case. 

There should be a change in the haze profile near 20 km, 
where the vertical optical depth of methane is equal to 0.5 at 
1216 A. At the slant optical depth t CH 4 > 1 below this level, 
the production of the condensible species becomes low, and 
the haze and its specific surface are distributed with the 
atmospheric scale height H due to variation of the particle 
sedimentation velocity. Methane becomes optically thin 
above 20 km, and the haze specific surface is distributed 
with a scale height 

H<m ») 

The mean methane scale height at the two observed occulta- 
tion points is 7 km, and H s = 7 km as well. We have found 
a function which approximates rather well the condensation 
rate coefficient above and below 20 km. 

The condensation rate coefficient adopted by LYA93 at 9 
km is equal to 10~ 5 s’ 1 and is smaller by a factor of 4.7 than 
our value, which is based on the measurements. Fortunately, 
they assumed a haze scale height of 10 km, which is just the 
mean of our values below and above 20 km. 

Boundary conditions. Photochemical processes are 
strongly depleted in the troposphere (below 8 km) due to 
high opacity of the overlying atmosphere to the solar 
radiation at A<1400 A, and due to condensation with a time 
constant which is much smaller than the mixing time. 
Therefore condensation of photochemical species on the 
cloud particles (which exist mostly below 8 km [ Rages and 
Pollack , 1992]) and on Triton’s surface is of minor impor- 
tance and is neglected in our model. 

We assume in our basic model that the haze particles are 
chemically passive. An extreme alternative case, when 
scavenging of all species except N 2 , CH 4 , CO, and H 2 occurs 
with efficiency of 100%, will be also considered. 

H, H 2 , C, N, and O are subject to thermal escape at the 
upper boundary at 800 km. This boundary is slightly lower 
than the exobase height of 870 km [Krasnopolsky et al., 
1993]; however, this fact is of minor importance to the 
results of modeling. All other neutral compounds are 
assumed to have zero fluxes at the upper boundary. 

SS95 argue that Triton’s ionopause is rather close to its 
exobase. Ions reaching the ionopause are swept out by the 
corotating magnetosphere of Neptune. This process is 
responsible for the observed decrease of the plasma scale 
height above 600 km. To account for this process, Yung and 
Lyons [1990] and Summers and Strobel [1991] assumed an 
ion escape velocity v- * 5*10 4 cm/s, i.e., "at the maximum 
ambipolar diffusion rate. 11 This velocity is too high and 
would severely deplete the C + ionosphere. Therefore LY A93 
adopted a minimum ion escape velocity, which is the thermal 
escape velocity of 10 cm/s for C + . Though small v i facilitates 
modeling of the electron peak density, we assume v, = 150 
cm/s in our basic model. This velocity reflects the possible 
ion escape above the ionopause and helps to reproduce the 
electron density profile above 600 km. 

The diffusion coefficients of neutral species in our models 
are based on the data taken from Chamberlain and Hunten 
[1987] and Banks and Kockarts [1973]: I\ = 4.87* 10 17 
7 €698 /«, D m = 1.88* 10 l7 7* 82 //t, D N c 0 = 10 ,7 7 ° 75 /*, and D = 
0.7* 10 17 7® 75 //i for other neutral molecules. Ambipolar 
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diffusion coefficients of ions are calculated using a method 
described by Atreya [1986]: Z^ + = 5>< 10 18 //?, D c+ N+ 0+ = 
10 i8 //j. All these values are in cm 2 /s, and n is the total 
number density. 

Continuity equations in a spherical atmosphere are solved 
in our model by a method described by Krasnopolsky 
[1993a] and Krasnopolsky et ai [1993]. The accuracy of the 
finite difference method we use depends strongly on the 
number of steps. The simplest check is the application of this 
method to species which do not participate in chemical 
processes. In this case, a decrease in the number density up 
to height h is given as (1 - hfHf® instead of Q h/H . Here H is 
the scale height, assumed constant for the sake of simplicity, 
and 6 is the altitude step. Then a ratio of the calculated 
density at the upper boundary to the real one is equal to 
0.01! for 50 total steps, 0.145 for 100 steps, 0.29 for 200 
steps, and 0.58 for 320 steps. We feel that 58 steps used by 
LYA93 do not provide sufficient accuracy with our method, 
and 320 steps will be used. 

5. Eddy Diffusion 

Methane profiles from the UVS solar occultation measure- 
ments [Herbert and Sandel, 1991] are used to obtain the 
eddy diffusion coefficient K. There are two competing 
processes which determine the CH 4 profile below 100 km on 
Triton: upward transport by eddy and molecular diffusion 
and photolysis by solar 1000-1400 A and Ly man-alpha 
background radiation. The idea of explaining the measured 
CH 4 profiles by these processes was suggested by Strobel et 
ai [1990a], who applied it to preliminary data for CH 4 and 
obtained K - 6000 ± 2000 cm 2 /s for a quantum yield of CH 4 
photolysis y = 0.6. SS95 repeated this procedure with the 
improved methane profiles from Herbert and Sandel [1991] 
and gave K = K 0 (n 0 ln) l/2 cm 2 /s, where K^ = 1600 and 1200 
cm 2 /s at the egress and ingress occultation points, respective- 
ly, and n is the atmospheric gas number density. We modi- 
fied their method to our model with y = 1.1. The balance of 
the methane is determined by its photolysis with column rate 
(CR) 2.25 x 10 8 cm'V 1 (Table 1), its reaction with CH ((R26) 
CR = 9.5* 10 7 cm'V 1 ), and by formation of CH 4 in ((R20) 
CR = 7.2* 10 7 cm'V 1 ). To account for these processes, an 
effective yield of methane photolysis is equal to 1.1. We 
assume K = 10 5 cm 2 /s in the troposphere from Yelle et ai 
[1991], though this does not improve agreement between the 
calculated and the measured methane profiles. The measured 
methane densities near the surface are uncertain within a 
factor of 4, and we will therefore use them as the second 
fitting parameter (the first one is K). Local values of the 
orbit-mean solar flux are used for these calculations. The 
values of K obtained are rather close for both occultation 
sites, and we therefore assume them to be equal. The best fit 
to the measurements is for either K - 2000 (n/nQ) I/2 cm 2 /s or 
K = 4000 cm 2 /s. The latter fits slightly better and will be 
used in our models (Figure 1). Then the homopause is at 32 
km, and methane profiles calculated for different K become 
parallel above the homopause (Figure 1). Our K are rather 
close to those of SS95, which located the homopause at 35 
km (40-50 km given in SS95 is not correct). The divergence 
of the CH 4 profiles calculated for different K is much 
smaller in Figure 1 than in SS95. Our calculations favor a 
CH 4 density of 3. lx 10 cm near the surface. 



Figure 1. Comparison of the measured profiles of methane 
[Herbert and Sandel , 1991] with the calculated profiles for 
[CHJ = 3.1xl0 n cm 1 near the surface and eddy diffusion 
coefficient K = 2000, 4000, and 8000 cm 2 /s. K = 10 5 cm 2 /s 
is assumed below 10 km [Yelle et ai , 1991]. The summer 
profiles are displaced to the left by a factor of 10. 


6. Basic Model: Neutral Atmosphere 

Photon and electron dissociation and ionization. 

Among the reactions in Tables 1 and 2 are 13 dissociation, 
10 ionization, and four dissociative ionization processes. To 
calculate these processes we use absorption and ionization 
cross sections of N 2 and O from Kirby et ai [1979], of CO 
(for X<800 A) from Torr et ai [1979], of N from Comes and 
Elzer [1968], and of C from Cantu et ai [1981]. Photolysis 
rates of CH 4 and other hydrocarbons and HCN may be 
calculated using their known cross sections and quantum 
fields (see Table 1). The most difficult ranges are 800-1000 
A, where absorption by N 2 bands dominates, and 1000-1 100 
A, where CO and C absorb. 

Compilations of the N 2 cross sections in the 1 A interval 
by Kirby et ai [1979] cannot help much at 800-1000 A 
because the line widths are smaller by 3 orders of magnitude 
than 1 A. Detailed modeling of the N 2 bands based on the 
best recent measurements with spectral resolution of 5 mA 
[Stark et ai , 1992] should be done with special care to solve 
the problem of absorption at 800-1000 A. An approximation 
which may be rather good for our purpose is based on the 
Voyager UVS occultations of Triton which gave the mean 
optical thickness of N 2 in this range as 

t = 1.1x10' 8 {N 2 )° 434 

where {N 2 } is the N 2 column density in cm' 2 [Krasnopolsky 
et ai , 1993]. The N 2 dissociation rate is equal to 


p(z) = 


2 dz 


where / is the total solar intensity in 800-1000 A at the orbit 
of Triton. The rate is maximum at 210 km, and the calculat- 
ed atmospheric transmission at 800-1000 A is equal to 0.1 at 
200 km for the geometry of the solar occultation. This 
transmission agrees with the measured transmission [Krasno- 
polsky et ai , 1993, Figure 1]. LYA93 applied the data from 
Kirby et ai [1979], finding that the calculated N 2 dissocia- 
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tion rate is maximum at 470 km. This corresponds to a 
completely opaque atmosphere for solar occultation at 350 
km and disagrees with the measured transmission of 0.5 at 
this height (see the same figure). 

Even more difficult is the calculation of the dissociation 
of CO at 885-1118 A. The bulk of the dissociation is in the 
region 1000-1118 A, which is not contaminated by N 2 
absorption. Four CO bands centered at 1003, 1052, 1062, and 
1076 A absorb in this range [Fox and Black , 1989]. We 
calculated a position of the maximum photolysis near 20 km. 
A strong methane absorption sharply reduces the CO 
dissociation below 20 km. The column dissociation rate is 
equal to 3.8 X 10 5 cm'V 1 , which is smaller by a factor of 3 
than that of LYA93. The CO photodissociation is a minor 
process in the CO destruction, and uncertainty of the 
photodissociation rate is not very important. 

Electron energy deposition in N 2 was considered by Fox 
and Victor [1988], and we use their data. For the sake of 
simplicity, we assume that the CO electron dissociation, 
ionization, and dissociative ionization cross sections are 
equal to those of N 2 , and N, O, and C electron ionization 
cross sections are equal to half those of N 2 . These processes 
are of little importance, and these approximations do not 
affect the results of modeling. All of the primary processes 
in the N 2 atmosphere are summarized in Table 3. Their 
vertical profiles are very similar to those calculated by SS95 
and are not given here. 

Hydrocarbons and HCN. The calculated profiles of the 
seven most abundant species are shown in Figure 2, the 
profiles of hydrocarbons and HCN are given in Figure 3, and 
some active radicals are shown in Figure 4. 

One half of the CH 4 photolysis events produces methyl 
CH 3 . Two thirds of these radicals recombine to CH 4 and the 
rest to ethane C 2 H 6 . These three-body reactions occur mostly 
below 50 km, and almost all ethane condenses on the haze 
particles with a precipitation rate of 28 g/cm 2 /b.y. The other 
half of the CH, photolysis events forms CH, which reacts 
with CH 4 (R26). The reaction product, ethylene C 2 H 4 , 
condenses with a rate of 135 g/cm 2 /b.y. Two percent of the 
C 2 H 4 photodissociates and reacts with C (R52) to form 
acetylene C 2 H 2 . Acetylene condenses mostly below 50 km 
with a rate of 1.3 g/cm 2 /b.y. and reacts with C (R51) above 
50 km. Both and C^ have similar production and loss 
rates above 50 km ((R52) and R(51)); therefore their profiles 
coincide in this region. A sixth of the CH radicals forms 
HCN via the (R24), (R60) cycle (see section 3), which starts 
from the reaction of CH with N 2 . Obviously this source of 
HCN is stronger (by a factor of 50 in this model) than the 


Table 3. Column Production Rates of N 2 \ N\ and N 


Process 

Column Rate 

Magnetospheric electron ionization' 


6.9 

Photo- and photoelectron ionization 


2.9 

Photo- and photoelectron dissociative ionization 

0.42 

Photodissociation 800-1000 A 


1.8 

Photoelectron dissociation 


0.91 

Production of N 

15 

(5.8) 

Production of N 2 + 

8.6 

(2.9) 

Production of N* 

r» * irt/ -Z 1 

1.7 

(0-42) 

.L _ C 


Rates are in units of 10' cm'V 1 normalized to the surface. 
Numbers in parentheses are for the solar model. 


* Production of N 2 + 82%, production of N + 18%, electron impact 
dissociation 60% of this value. 


20 40 60 60 100 T (K) 



Figure 2. Composition of Triton’s atmosphere: the most 
abundant species. The N 2 and T profiles and the measured N 
profile (short dashes), N 2 , and T at 575 km are from Krasno- 
polsky et al [1993]. 

reactions N + CH 3 , CH 2 ((R27), (R30), and (R59)). HCN is 
rather inert, and its sinks are condensation below 50 km with 
a rate of 29 g/cm 2 /b.y. and reactions with ions above 200 
km. Coincidence of the HCN and QH* profiles below 40 km 
(Figure 3) is accidental. A photosensitized dissociation of 
CH 4 [Yung et al . , 1984] via 
(R2) C 2 H 2 + /iv - C 2 H + H 
(R70) C 2 H + CH 4 - C 2 H 2 + CH 3 
net CH 4 + /iv - CH 3 + H 

is strongly depleted due to condensation of C 2 H 2 and loss of 
C 2 H in the reactions with N and O. 

Atomic and molecular hydrogen. According to section 
2, the production of hydrogen in both forms is equal to 
2.25 x 10 8 x2x 1.1 = 5* atoms/(cm 2 s). Here the first term 
is the total global -mean 1000-1400 A photon flux, two 
hydrogen atoms are released per photon, and 1.1 is the yield 
of methane photolysis. Hydrogen has no important sinks 
except thermal escape; therefore we know a priori the total 
hydrogen escape rate. Atmospheric and ionospheric reactions 
result in mostly mutual transformations between atomic and 



Figure 3. Composition of Triton’s atmosphere: photochemi- 
cally produced hydrocarbons, HCN, and CH 3 . 
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Figure 4. Composition of Triton's atmosphere: diatomic and 
triatomic radicals. 


Table 5. Main Production and Loss of H 2 


Process 

Column Rate 

CH 4 + hv (Y = 0.514) 

1.32 

h + h + m-»h 2 + m 

0.68 

NH + H - N + H 2 

0.46 

Total production 

2.46 

Ion reactions 

0.46 

Escape 

1.96 

Total loss 

w-.: • . -C* 1 rX TT- 

2.42 

l 


Sources and sinks of N are given in Table 6. The structure 
of the profile below 50 km is due to indirect photolysis of N 2 
via the CH+N 2 cycle discussed in section 2. Atomic nitrogen 
loss is due to the three-body association reaction (R1 1) and 
catalytic cycles with C, O, and H. The C cycle, which 
includes the reactions (R25), (R41), and (R45), was dis- 
cussed in section 3. The O cycle is represented by the 


molecular hydrogen with a small net effect on the sum of 
these species. Molecular hydrogen is formed from H atoms 
below 100 km by the three-body reaction (R12) and by the 
reaction (R37) NH + H -*■ N + H 2 . NH is produced by the 
CH+N 2 cycle and by the three-body reaction (R16). Reac- 
tions of N 2 + , N + , and CO + with H 2 result in its indirect 
photolysis. Balances of production and loss of atomic and 
molecular hydrogen (normalized to the surface) are shown in 
Tables 4 and 5. These tables as well as similar tables for 
other species discussed below include only the main process- 
es, so the balances of production and loss are not exact. 
Complete data for the column reaction rates from Tables 1 
and 2 may be used in more careful calculations if needed. 
The total hydrogen escape coincides with its production. The 
H 2 mixing ratio is rather constant with height due to the 
diffusion-limiting flow. The H profile is maximum at 25 km. 
A factor of 5 decrease at the surface shows that both 
downward flux of H and locally produced H atoms recom- 
bine almost totally below 50 km. 

Atomic nitrogen* The calculated profile of N is com- 
pared with the measured profile [Krasnopolsky et al, 1993] 
in Figure 2. The measured profile has number densities (1 ± 
0.25)* 10 8 cm' 3 at 400 km and (5 ± 2.5) * 10 8 cm' 3 at 200 km, 
which agree with the calculated densities 7.2* 10 7 cm 3 and 
5.7* 1QP cm' 3 , respectively. Krasnopolsky et al. claimed that 
the profile should correspond to diffusion equilibrium above 
470 km, where the measurements were rather uncertain. 
SS95 are doubtful of this extrapolation. The calculated 
profile shows that the extrapolation from 470 km to the 


reactions 

(R55) N + O NO + hv 
(R57) N + NO -► N 2 + O 
net N + N — ► N 2 

The H cycle is similar, except for the three-body formation 
of NH. The oxygen and carbon catalytic cycles are of 
comparable importance, and each provides one third of the 
atomic nitrogen loss each. The carbon cycle was dominant in 
LYA93 (96%) and minor in SS95. The smaller role of the 
carbon cycle in our model than in LYA93 may be due to a 
weaker temperature dependence adopted for and smaller 
k 4l in our model. Another loss mechanism for N is thermal 
escape, which is equal to 7.7* 10 24 s' 1 and agrees with the 
measured value of (1 ± 0.3)* 10 25 s' 1 . 

Carbon monoxide, atomic oxygen and carbon. As 
discussed above, the carbon monoxide loss rate by direct 
photodissociation and ionization is 3.8 x 10 5 cm V. A charge- 
exchange reaction of CO with N 2 + is of comparable impor- 
tance (4.8* 10 5 cmV 1 ), while that with 1ST is the dominant 
loss process with a column loss rate of 1.5* 10 7 cm~ 2 s At a 
rate of 6*10 6 cmV 1 , CO molecules are given back to the 
atmosphere by charge exchange of CO + with H 2 , C, O, H, 
CH 4 and by recombination of HCCf. Carbon monoxide is 
formed by reactions of atomic oxygen with CNN and CN 
(4*10 6 cm'V 1 ) and with CH, radicals (5*10 6 cm'V 1 ). The 
CO mixing ratio in the upper atmosphere is smaller than near 
the surface due to large ionospheric losses. The total CO loss 


exobase is accurate within 5%. 


Table 6. Main Production and Loss of N 



Process 


Column Rate 

Table 4. Mam Production and Loss of H 

N, + (/tv, e) 

- N, NT 

1.70 

Process 

Column Rate 

The CH + R, cycle 

0.20 

CH 4 + hv (y = 0.996) 

2.24 

(N 2 + + (e, 0»x2 

0.87 

CH 4 + CH - + H 

0.94 

Total production 

2.77 

Ion reactions 

0.93 




Total production 

4.11 

(N + NO -► 

N 2 + 0)*2 

LOO 

(H + H + M-+H 2 + M)x2 

1.37 

(N + CN -► 

N 2 + C)*2 

0.90 

H + CH 3 + M — CH 4 + M 

0.72 

(N + N + M 

-» N 2 + M)x2 

0.33 

(NH + H - N + H)*2 

0.92 

(N + NH -► 

N 2 + H)*2 

0.21 

Escape 

1.04 

Escape 


0.34 

Total loss 

4.05 

Total loss 


2.77 


Values are in units of 10 8 cmV 1 . Values are in units of lCf cmV. 
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is balanced by its production with a small difference due to 
the oxygen escape at the rate of 1.9* 10 5 cmV, and due to 
formation and precipitation of C0 2 at 5.3* 10 3 cmV. This 
difference is supplied by evaporation of the surface ice. The 
net effect of the photochemistry is the loss of 0.3 g/(cm 2 b.y.) 
of carbon monoxide. The production and condensation of 
C0 2 are extremely low and cannot be responsible for the 
observed C0 2 in the ice [ Cruikshank et al , 1993], 

Atomic oxygen (Table 7) is formed by dissociative 
recombination of CO + and by photolysis of CO mostly at 
1000-1080 A. The latter has a strong cutoff below 20 km 
where methane absorbs. We assumed no CO photolysis 
below 20 km, and the profile structure in the lowest 50 km 
reflects the effects of the CO photolysis. Oxygen is removed 
mostly in the reactions with CFI* and CN X radicals. 

The sources of atomic carbon (Table 8) are CH 4 , CO, and 
HCN photolysis (the latter ends by (R45) CN + N), reactions 
of CH with N and H ((R32), (R33)), and recombinations of 
CO + and some other ions. Another source is the cycle 
(R51) C 2 H 2 + C C 3 + H 2 
(R65) C 3 + N - CN + C 2 

(R63) C 2 + N -► CN + C 
(R45) (CN + N -> N, + C)x2 

net C 2 H 2 + 4N - 2N 2 + 2C + H 2 
This cycle is questionable because products of (R51) are not 
known, and those suggested produce the highest exothermi- 
city. This cycle provides 10% total production of carbon 
atoms and is not of crucial importance for the model. Carbon 
is lost in reactions of CNN with O and H (because CNN is 
formed by C + N 2 + M and mostly gives carbon back in 
CNN + N and CN + N), in (R53) C 2 H 4 + C, and in charge 
exchange with N 2 + and CO + , with minor contributions from 
some other processes. A quarter of the C production escapes. 
Escape of carbon is greater by a factor of 25 than that of 
oxygen (due to the difference in their masses) and is supplied 
by methane. 

Diatomic and triatomic radicals. The most important of 
these radicals are shown in Figure 4. Many of these species 
have maxima at 60 km caused by the reactions (R5I), (R52) 
C 2 H 2 , C 2 H 4 + C. Carbon atoms (and CH 2 ) are released rather 
than lost in these reactions; this fact impacts the densities of 
all carbon-bearing radicals and is seen even in the C profile 
(Figure 2). 

7. Ionosphere 

The column ionization rate normalized to the surface is 
equal to 10 8 cm V 1 (Table 3), with a ratio of solar to 


Table 7. Main Production and Loss of O 

Process 

Column Rate 

ccr + -* c + o " 

9.19 

CO + hv -» c + o 

0.38 

Total production 

9.57 

O + (CH, CH 2 , CH 3 , C 2 H 4 ) 

5.12 

O + (CN, CNN) 

3.98 

O + (Q, C 3 ) 

0.16 

O + O + M— >0 2 + M 

0.15 

Escape 

0.19 

Total loss 

9.60 

Values are in units of 10° cm' 

V. 


Table 8. Main Production and Loss of C 

Process 

Column Rate 

CO + hv — C + O 

0.38 

CH 4 + hv -» C + H 2 + H 2 

0.90 

CH + H - H 2 + C 

2.82 

CH + N -+ H + CN’ 

1.91 

(CjHj + C)«2 

2.04 

CO + + e -» C + O 

9.19 

Other ion recombinations 

0.61 

Total production 

17.8 

C + C 2 H 4 - C 2 H 2 + ch 2 

1.47 

CNN + + (H, O) - N 2 + (CH, CO) 

9.00 

C + (N 2 \ CO + ) -» C + + (N 2 , CO) 

1.42 

Escape 

4.68 

Total loss 

16.6 


Values are in units of 10° cm'V 1 . 

* The reaction CN + N -► N 2 + C ends this cycle. 
f These cycles begin from C + N 2 + M — CNN + M. 


magnetospheric electron rates of 1:2. If the observed iono- 
sphere, with a peak electron density of (3.5 ± 1)* 10 4 cm 3 at 
340 km, consisted of molecular ions, it would require an 
ionization rate larger by 2 orders of magnitude than the 
above value due to rapid dissociative recombination of 
molecular ions. Therefore atomic ions having very slow 
radiative recombination rates dominate the ionosphere. 

Lyons et al [1992] established the importance of atomic 
carbon in the formation of the ionosphere. C + is formed by 
charge exchange between N 2 + and C (R84) and is removed 
by reactions with CH 4 and HCN ((R105), (R106)). C + ions 
dominate in their model of the ionosphere, their lifetime is 
rather long, and the calculated electron densities fit the 
measured values even without ionization by magnetospheric 
electrons. However, the reaction (R84) rate coefficient must 
be very large {k% A = 10' 9 cm 3 /s) to facilitate production of C + 
ions. 

The introduction of CO has many important aspects to the 
ionosphere. The reaction with CO (R90) becomes a principal 
sink of N + , substantially lowering its densities. Recombina- 
tion of CO + increases the atomic carbon densities at 200-400 
km, and it becomes possible to fit the measured electron 
densities with reasonable values of k 84 95 . CO\ HCO\ and 
NO + become important molecular ions. FT ions are supressed 
by rapid charge exchange with atomic oxygen (R99), with 
subsequent reaction between Q + and N 2 . Reaction (R99) is 
slightly endothermic, and we therefore reduce its rate 
coefficient to a value appropriate at 100 K. The time 
constant for C + ions above 250 km is much longer than a 
period between passes of the magnetic equator (7 hours), and 
temporal variations of electron density become small. This 
corresponds to the global mean and orbit-averaged conditions 
typical of our model. 

In our model, C + is a dominant ion (Figure 5), and its 
balance is considered in Table 9. As in Lyons et al [1992], 
the main source of C + is charge exchange between N 2 + and 
C. The rate coefficient of this reaction is a free parameter to 
fit the measured electron density profiles. The contribution 
of the similar reaction with CCf instead of N 2 + with the same 
rate coefficient is smaller by an order of magnitude. The 
total production of C + is only about 2% of the total ioniza- 
tion rate. The chemical loss is by reactions (R105), (R106) 
with CH 4 and HCN, and radiative recombination. The 
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strongest sink in our model is the ion escape. As discussed 
in section 4, the escape velocity is unknown and is chosen as 
150 cm/s to better fit the shape of the electron density profile 
between 400 and 700 km. Besides, if the ionopause is really 
close to the exobase (SS95), then the ion escape velocity 
should substantially exceed the thermal escape velocity. If 
the assumed escape velocity is much smaller than 150 cm/s, 
then the calculated electron density at 700 km becomes 
larger by a factor of 2. A significant reduction of the ion 
escape velocity makes it possible to reduce k^ 95 by a factor 
of 3 with the same electron densities below 400 km, while 
the densities near 700 km become larger by a factor of 2. 

In our model, hT and C + densities are comparable above 
500 km, where loss of N* is small due to low densities of 
CO and H*, and are similar to SS95 in this respect. While 
the bulk ionosphere is of an F type, the region below 240 
km consists of molecular ions and is thus of E type. As 
discussed above, the primary ions N 2 and N react mostly 
with either H 2 forming N 2 H + (in the case of N* , NH* ions 
formed react with N 2 to produce N 2 H + ) or CO to yield CCf . 
The main sinks of N 2 H + and CO + in the lower ionosphere are 
reactions with CO and H 2 , respectively. A product of these 
reactions, HCO\ becomes the major ion in the E region. 

8, Other Models 

Some parameters of our model are rather uncertain, and 
we have therefore calculated models with different values of 
the parameters. The main properties of these models are 
shown in Table 10. Model 1 is the basic model without the 
CH + N 2 cycle (k m = 0 for this model). Evidently, formation 
and sedimentation of HCN becomes smaller by a factor of 50 
than that in the basic model. More CH radicals react with 
CH 4 to form QH 4 , and H atoms released in this reaction 
react with CH 3 reducing the production of ethane. The other 
principal properties of the model remain unchanged. 

Model 2 is our basic model without the reaction (R25) C 
+ N 2 + M. The formation of CNN affects loss of C and O 
via reactions (R42), (R44) with H and O. Carbon and oxygen 
densities are larger in model 2 than in the basic model. This 
results in a small increase in the C and electron densities. 
The loss of nitrogen becomes smaller (due to the absence of 


the CNN cycle) and larger (due to increase of the oxygen 
cycle; see Table 6) with a net increase by only 4%. 

Chemically passive haze particles were assumed in the 
basic model, and condensation of only C 2 H X and HCN was 
allowed to occur on these particles. Model 3 assumes that all 
species except N 2 , CH 4 , CO, and H 2 are scavenged by the 
haze at the maximum possible rate. The strongest effect of 
this assumption is an increase of ethane production by a 
factor of 3 due to scavenging of H atoms and smaller loss of 
CH 3 in (R20) CH 3 + H. Hydrogen atoms do not form H 2 , 
reducing the H 2 densities and total hydrogen escape flux by 
a factor of 2. Reaction (R82) N 2 + + H 2 is the main loss 
mechanism for N 2 + ; therefore N 2 + and production of C + are 
larger due to the reduction of H 2 in this model. The peak 
electron density is larger than the measured value by a factor 
of 1.5, and a reduction in k gAt 95 may be applied to fit the 
measured electron density. 

The main properties of our model are rather stable even 
under the extreme conditions of models 1, 2, and 3. Now we 
consider variations of the model results arising from changes 
in the assumed CO mixing ratio. These changes are compen- 
sated by a proper choice of k% A 95 to fit the measured peak 
electron density. Model 4 has f co = 10 4 and A: 84 95 = 8 X 10 
cm 3 /s (Figure 6), and model 5 is for 10 3 and 2 x lCf 11 cm/s, 
respectively. The effects of the CO variation (by a factor of 
10) are nonlinear: variations of the CO density at 300 km are 
a factor of 24, and variations of C and O densities above 250 
km and of their escape fluxes are a factor of 5. The product 
fcoksA, 95 ts not constant. 

Comparing model 4 with our basic model, we see the 
expected reductions in the densities of CO, C, and O 
(neutrals) and of CO + and NO + (ions). N , N 2 , N 2 FF, and H 
densities are rather insensitive to the CO change. 

Model 6 is our model adjusted to the conditions of the 
model of LYA93: solar (without magnetospheric electrons) 
ionization, the CO mixing ratio is 10 4 , the ion escape is 
thermal (v, = 10 cm/s), the charge-exchange rate coefficient 
is 10' 10 cm 3 /s, the CH + N 2 cycle is absent, and the radiative 
recombination coefficients are 5><10 12 cmVs. Despite differ- 
ences in other reactions and their rate coefficients and in the 
aerosol properties (see section 4) between model 6 and 
LYA93, the results are rather similar. This confirms the idea 
of the stability of the model relative to variations in input 
data, except for a few (the magnetospheric electron energy 
input ,/ c0 , k u> 95 > V i)* 

Model 7 is adjusted to that of SS95: the magnetospheric 
electron energy input is 1.4x10 s W, the CO mixing ratio is 
2.5 x lO 4 , k* 4 95 = 10 “ cra’/s, v, = 10 cm/s (we do not know 
their v ( ). Again the electron density profile in model 7 is 

Table 9. Main Production and Loss of C* 


Process ^ Column Rate 

N-f + C -* C + N, L64 

CO* + C-C* + CO 0.19 

C + hv C* + e 0.07 

Total production L90 

C* + e — C + /iv 0.23 

C + CH 4 -* products 0.37 

C* + HCN -» C.N* + H 0.17 

Escape L14 

Total loss 1-9 1 


Values are in units of 10° cm V. 
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Table 10, Models of Triton’s Atmosphere and Ionosphere 


Value 

Measured 

Basic 

1 * 

2 T 

3 * 

4 

5 

“6* 

T 

fee *10 ’ 

- 

3 

3 

3 

3 

1 

10 

1 

2.5 

*(rv,co‘+Q, *io 11 

- 

4 

4 

4 

4 

8 

2 

10 

1 

v,, cm/s 

- 

150 

150 

150 

150 

150 

150 

10 

10 

/* fnM , km 

340 

320 

320 

300 

300 

290 

330 

280 

290 

x 10 4 cm 5 

3.5 ± 1 

3.5 

3.6 

4.25 

5.1 

3.6 

3.0 

3.3 

3.4 

*700 km. *10 4 cm 

0.55 ±0.1 

0.76 

0.75 

0.78 

1.01 

0.77 

0.56 

1.00 

1.52 

I N *l 700 km> *10" cm' 3 

- 

0.34 

0.34 

0.33 

0.57 

0.45 

0.13 

0.09 

0.71 

[Nhoo km, * 10 8 cm 3 

5 ±2.5 

5.7 

5.7 

5.9 

5.9 

6.5 

4.9 

2.5 

7.5 

|NW„> *10 7 cm 3 

10 ±3 

7.2 

7.2 

7.5 

7.5 

7.7 

6.5 

3.1 

9.6 

ICOImow *10^ cm 3 

- 

1.1 

1.1 

1.0 

0.7 

0.3 

6.5 

0.6 

0.7 

ICJ300 lon> * cm 3 

- 

3.9 

3.8 

4.3 

4.4 

1.4 

7.1 

0.9 

4.0 

<I> N> *10 24 s' 1 

10 ±3 

7.7 

7.7 

8.0 

8.2 

8.6 

6.8 

3.2 

10.0 

«I>H2. X 10^ s’ 

- 

4.5 

4.6 

4.8 

2.1 

4.4 

4.8 

5.4 

4.2 

4> h , * 10 25 s' 1 

- 

2.4 

2.4 

2.2 

1.6 

2.7 

1.8 

0.7 

3.3 

<J> 0 , *10“ s' 1 

- 

4.4 

4.4 

5.1 

4.8 

1.7 

8.4 

0.8 

4.3 

<D C , * I0 24 s 1 

- 

1.1 

1.1 

1.2 

1.2 

0.4 

2.0 

0.24 

1.1 

<t> ( „ * 10 23 s' 1 

- 

2.6 

2.6 

2.8 

2.8 

2.0 

2.7 

0.4 

0.4 

S nic , g cm 2 b.y.' 1 

- 

1.3 

1.5 

3.3 

0.7 

1.2 

0.8 

1.3 

1.4 

g cm' 2 b.y.' 1 

- 

135 

164 

133 

155 

137 

134 

166 

166 

S’.ciio. g cm' 2 b.y.' 1 

- 

28 

22 

28 

81 

29 

27 

22 

22 

•S’hcti. g cm 2 b.y. 1 

- 

29 

0.57 

36 

10 

32 

28 

3 

3 


Hereto i g the CO mixing ratio, v i is the ion escape velocity, h trvut and e ^ are the altitude and electron density at the 
ionospheric peak, <f> ( and S t are the total escape and sedimentation fluxes of i species. 

* Basic model without (R60) HCN 2 + H. 

f Basic model without (R25) C + N 2 + M. 

* Basic model with condensation and scavenging of all species except N 2 , CH 4 , CO, and H 2 . 

5 Analog of the solar model of LYA: (R60) HCN 2 + H is neglected, radiative recombination coefficients are 5* 10 12 
cm 3 /s. 

1 Analog to SS with P = 1.4*10 8 W and without (R60). 


rather similar to that of SS95 (we do not have other results 
of SS95). 

Of all our models, only model 6 fails to reproduce the 
measured profile of atomic nitrogen. This is the solar model 
(without precipitation of magnetospheric electrons). There- 
fore the main conclusions of the model comparison are that 
(1) magnetospheric electron energy input of about 10 8 W, 
i.e., approximately twice the solar EUV input, is needed to 
reproduce the measured N densities, and (2) it is impossible 
to establish firm constraints to the CO mixing ratio by 
modeling the atmosphere and the ionosphere with the 
unknown rate coefficient /^ 4 . The plausible range 10 1 1 - 10' 10 


cm Vs for k% A requires the CO mixing ratio in the range 10 4 - 
10 \ Measurements of the CO mixing ratio in Triton's atmo- 
sphere, and laboratory measurements of the charge exchange 
between N 2 + and C are of great importance to Triton's 
photochemistry. 

9. Conclusions 

We have considered the photochemistry of 32 neutral and 
21 ion species in Triton's atmosphere. The model shows that 
three parent species N 2 , CH 4 , and CO (with a mixing ratio of 
3 X 10 4 in our basic model) sublime from the ice with rates of 




Figure 6. Model 4 with the CO mixing ratio of 10* 4 and k( N 2 + + C) = 8* 10 11 cmVs: basic (a) neutral and 
(b) ion species. 
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40, 208, and 0.3 g/cm 2 /b.y., respectively. Photolysis of 
methane by the solar and the interstellar background Lyman- 
alpha radiation drives the chemistry below 50 km and 
produces hydrocarbons C 2 H 4 , C 2 H*, and C 2 H 2 . These hydro- 
carbons form the haze particles and precipitate with rates of 
135, 28, and 1.3 g/cm 2 /b.y., respectively. The CH + N 2 cycle 
we introduced increases the production of HCN by more than 
an order of magnitude to a value of 29 g/cm 2 /b.y., and results 
in the indirect photolysis of N 2 . Using the method proposed 
by Strobel et al. [1990a], we fmd from the measured 
methane profiles the eddy diffusion coefficient K = 4*10* 
cm 2 /s above the tropopause and a methane number density of 
(3.1 ± 0.8)* 10 11 cm' 3 near the surface. Solar EUV radiation 
(X < 1000 A) and precipitation of magnetospheric electrons 
are responsible for chemistry above 200 km. We assume a 
magnetospheric electron energy input of 10 8 W based on the 
thermal balance calculations. N, H 2 , H, O, and C are the 
most abundant photochemical species, with total escape rates 
of 7.7* 10 24 s' 1 , 4.5 x 10 25 s' 1 , 2.4* 10 25 s 1 , 4.4* 10 22 s' 1 , and 
1.1* 1 0 24 s' 1 , respectively. The chemistry of a very interesting 
intermediate region at 50-200 km is driven mostly by the 
atomic species, which are transported to this region from the 
atmosphere above and below it. The ionosphere consists of 
the E region at 150-240 km and of the F region above 240 
km. A major ion of the E region is HCO\ C + dominates in 
the F region and forms the peak at 320 km. The ionosphere 
above 500 km includes almost equal densities of C + and FT 
ions. The calculated profiles of atomic nitrogen and electron 
densities agree with the measurements. 

Some other models with some reactions removed or 
changed, with changed properties of the haze particles (from 
chemically passive to active at a maximum rate), etc., were 
developed. These models show that some of these variations 
are relatively unimportant for the basic structure and compo- 
sition of the atmosphere and provide some local changes in 
some species. However, there are four basic unknown values 
which have strong impacts on the composition of the 
atmosphere and ionosphere. These values and their plausible 
ranges are the CO mixing ratio f co - 10' 4 -10’ 3 , the magneto- 
spheric electron energy input (1 ± 0.5)* 10 8 W, the rate 
coefficient of charge-exchange reaction N 2 + + C k = 10 n - 
10~ 10 cm 3 /s, and the ion escape velocity v,. » 150 cm/s. 
Various combinations of these values allow us to construct 
models which fit the measured N and electron densities. 
Therefore measurement of CO in Triton’s atmosphere (which 
we hope can be made by ground-based spectroscopy) and 
laboratory study of the N 2 + + C reaction are very important. 
Modeling shows that the measured profile of atomic nitrogen 
cannot be reproduced without a substantial magnetospheric 
electron energy input, which should exceed the solar EUV 
input by a factor of 2-3. If it does, then the ionopause is 
rather close to the exobase (SS95), and the ion escape 
velocity is much larger than that for thermal escape. This 
increases requirements on the N 2 + + C reaction. 

Some results discussed here may be used for modeling of 
the N 2 -CH 4 atmospheres of Titan and Pluto while others may 
not. For example, formation of HCN by the CH-N 2 mecha- 
nism is important only in atmospheres with a very low CH 4 
mole fraction, and the yield of HCN in methane photolysis 
may be given as 



where = 1.8*10 17 cm 2 is the cross section at 1216 A. 
This yield is negligible ( 10 8 ) on Titan, where f^ « 0.08 at 
the level of methane photolysis. 

Photochemical production of C0 2 is too low to explain the 
observed ratio of 10 3 in the ice [Cruikshank et al. , 1993], 
while the production of C 2 H X and HCN is too large to be 
missed in the observations. A reasonable explanation is 
decomposition of these species in the ice by the solar 
radiation X > 1400 A [Allen et al, 1994]. The atmosphere is 
transparent for this radiation, and other species may form, 
resulting in coloration of the ice [ Thompson and Sagan , 
1990]. Bohn et al. [1994] have suggested that formation of 
the highly reactive molecule CH 2 N 2 may synthesize larger 
alkanes from CH 4 and C 2 H 6 , driving organics in the surface 
ices toward saturation (i.e., to reduce multiple C-C bonds). 
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